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Abstract 

A major issue in general relativity, from its earliest days to the present, 
is how to extract physical information from any solution or class of so- 
lutions to the Einstein equations. Though certain information can be 
obtained for arbitrary solutions, e.g., via geodesic deviation, in general, 
because of the coordinate freedom, it is often hard or impossible to do. 
Most of the time information is found from special conditions, e.g., de- 
generate principle null vectors, weak fields close to Minkowski space (us- 
ing coordinates close to Minkowski coordinates) or from solutions that 
have symmetries or approximate symmetries. In the present work we 
will be concerned with asymptotically flat space times where the approx- 
imate symmetry is the Bondi-Metzner-Sachs (BMS) group. For these 
spaces the Bondi four-momentum vector and its evolution, found from 
the Weyl tensor at infinity, describes the total energy-momentum of the 
interior source and the energy-momentum radiated. By generalizing the 
structures (shear-free null geodesic congruences) associated with the al- 
gebraically special metrics to asymptotically shear-free null geodesic con- 
gruences, which are available in all asymptotically flat space-times, we 
give kinematic meaning to the Bondi four-momentum. In other words 
we describe the Bondi vector and its evolution in terms of a center of 
mass position vector, its velocity and a spin-vector, all having clear geo- 
metric meaning. Among other items, from dynamic arguments, we define 
a unique (at our level of approximation) total angular momentum and 
extract its evolution equation in the form of a conservation law with an 
angular momentum flux. 
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1 Introduction 



From the very earliest days of general relativity (GR) the issue of extracting 
the physical meaning or content in the solutions of the Einstein equations has 
been a difficult problem^l and usually has been solved only in special cases or 
special situations, e.g., in the presence of symmetries or approximate (asymp- 
totic) symmetries. Probably the best studied case is that of asymptotically flat 
space-times where an asymptotic symmetry exists, namely the BMS group, and 
where, at null infinity, the total energy-momentum of the interior source was 
identified as well as gravitational radiation i.e., energy and momentum loss. 
These identifications, in terms of the asymptotic components of the Weyl ten- 
sor, were made via group theoretical arguments combined with physical insight 
into the dynamics of the radiation process. However the problem of gravita- 
tional multipole moments, and in particular, the mass dipole and the angular 
momentum has proved to be difficult. At the present time, using symmetry 
arguments, there are several different proposed definitions [3J 01 [S]of these 
quantities with little apparent use of dynamic considerations. (To our knowl- 
edge only in linear theory or in stationary cases has the issue of quadrupole 
moments been addressed.) In the present work, using largely a dynamic ar- 
gument and verifying it by a symmetry consideration, we return to the issue 
of recognizing the mass-dipole moment and angular momentum hidden in the 
asymptotic Weyl tensor. From this dynamic argument the mass-dipole moment 
and angular momentum, using the Bondi energy-momentum four-vector, are 
found. Both kinematic expressions, as well dynamic equations (equations of 
motion), are determined. 

In flat-space Maxwell theory the related issues are much simpler though 
some questions do remain. The charge, obtained by a 2-D surface integral of 
the asymptotic Maxwell field is a constant and lies in the scalar representation 
of the Lorentz group; the dipole moments (electric and magnetic) also obtained 
by 2-D surface integrals of the asymptotic Maxwell field (though the results 
depend on a time direction and the 2-surface chosen for the integral) also lie 
in a finite dimensional representation of the Lorentz group. For a static or 
stationary system it is easy to define the center of charge (so that the electric 
dipole moment associated with it vanishes) and a constant magnetic dipole 
moment. Recently it was shown[6] that for general asymptotically flat Maxwell 
fields (with non- vanishing total charge) one could find a unique complex world- 
line in complex Minkowski space, referred to as the complex center of charge, 
from which the dynamic electric and magnetic moments could be obtained. 

It is this construction of the complex center of charge for the Maxwell field 
that we generalize to GR. This leads, in a very unorthodox manner, to a complex 
center of mass. The real part is identified with the usual center of mass while 
the imaginary part is the specific spin-angular momentum, i.e., the spin per 
unit mass. The argument leading to this complex world line is based largely 
on dynamics and analogies with electrodynamics and special solutions of GR, 
e.g., the algebraically special metrics and in particular, the charged spinning 
metric [7]. It is applied to both the asymptotic vacuum Einstein equations and 
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the Einstein-Maxwell equations. For ease of presentation most of the analysis 
will be for the vacuum case. The results for the Einstein-Maxwell case will be 
presented towards the end of Sec. Ill without detailed derivation. 

In order to get a better perspective on this work and on our point of view, 
several comments might be of some use. 

Both Maxwell theory and GR are considered to be fundamental physical 
theories [in addition to quantum theory which so far appears incompatible with 
GR] and thus in principle this construction is applicable to the gravitational 
field in the neighborhood of any massive body, from elementary particles to 
galaxies. Though GR is almost always applied to astrophysical situations, there 
is in principle no reason that it could not be applied at some level to laboratory 
masses or even to elementary particles. If the individual masses are sufficiently 
far apart we could even consider their gravitational fields at (relatively) large 
distances as being asymptotically flat and use the general theory of asymp- 
totically flat space-times for their analysis. The theory of asymptotically flat 
space-times appears to be the best tool to define isolated bodies in GR and, as 
in Newton's theory, the notion of isolated bodies is an approximation of reality 
that is useful in general, and specifically for defining center of mass, intrinsic 
angular momentum, etc. 

It is this point of view we adopt in this work. Though there is no rigorous 
way to make physical identification with the GR variables we will consider three 
very similar alternative identifications - some more intuitive than others - and 
see what are the theoretical consequences - and then argue for the last one. The 
idea is then to look at the dynamic predictions of each of the alternatives. At 
the linear level they all agree: it is in the consideration of the non-linear terms 
that we see the differences and see how one of the alternatives, probably the 
least intuitive one, comes closest to what we would physically expect. 

The starting point for the discussion are the asymptotic Bianchi identities. 
We have assumed that the Einstein or the Einstein/Maxwell equations have been 
integrated in the asymptotic region using the standard Bondi-Sachs-Newman- 
Penrose asymptotic [SJ HH [TUl [TT] (peeling theorem) behavior. When the relevant 
integrations have been completed we are left with several evolution equations, 
namely the asymptotic Bianchi Identities (the Bondi supplementary conditions) 
and asymptotic Maxwell equations [12] . that become our object for analysis and 
interpretation. 

In Sec. II we review the relevant ideas about asymptotically flat space-times 
(Einstein or Einstein/Maxwell) at null infinity[12] . More specifically we discuss 
the structure of null infinity (and its complexification) , i.e., Penrose's 3 + . We 
then review a variety of subjects: e.g., what 'lives' on 3 + , the asymptotic Weyl 
tensor with the Bianchi Identities and the Bondi shear, other related structures 
as cuts of 3 + , past-light-cones and angle fields. The spherical harmonic decom- 
position of the various fields is then described. In this context we discuss the 
Bondi four-momentum and the energy/momentum loss theorem. 

In Sec. Ill, we examine three alternative suggestions for extracting physical 
information from the Weyl tensor components in the asymptotic Bianchi Iden- 
tities. A central structure in this extraction process comes from the dynamics 
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associated with the Bondi four-momentum. More specifically we want first to 
identify the mass dipole (or center of mass) and angular momentum. Then, 
second, we want to establish dynamical laws for these quantities, i.e., we want 
equations of motion for the center of mass and evolution equations - with fluxes 
- for the angular momentum. They come directly from the energy-momentum 
loss equations and a certain reality condition. These dynamical laws are given 
at future null infinity and therefore are intended to apply to the entire isolated 
system. Interactions between different bodies can not be considered in this 
framework. 

The first and most 'obvious' or intuitive model for the selection for the mass 
dipole and angular momentum turns out to lack certain essential dynamic fea- 
tures and has little or no geometric meaning. We consider the first (and second) 
approaches as 'toy' models, introduced in order motivate and clarify the third 
and primary model. The second and the third approaches which are based on 
specific choices of complex slicing of future null infinity, referred to as (complex) 
cut functions, lead to very similar dynamics. The nod goes to the third choice, 
it being more 'natural' in its similarity or analogy to the algebraically special 
metrics and to its unity and geometric attractiveness. They are both based on 
a geometric construction of complex cut functions. One of them, the second, 
is more appropriate to flat space, the other (which does reduce to the second 
in flat space) is far better suited to the general asymptotically flat case. In 
the later case we give an extended description of the predictions and dynamic 
results including the results for the Einstein-Maxwell equations. These results 
include kinematic expressions for the Bondi mass and momentum, a center of 
mass and its equations of motion and a definition of angular momentum that 
includes both spin and orbital terms with a conservation law with an angular 
momentum flux. 

In Sec. IV we analyze the invariance, under the BMS group, of the third 
approach to this method of assigning physical meaning to the variables. We 
briefly discuss the representation theory of the Lorentz group due to Gelfand, 
Graev and Vilinkin [13l I14j . and in particular apply their description of the in- 
variant finite dimensional subspaces of the infinite dimensional representations, 
to our physical assignments. As a simple example, we apply the representation 
theory to the Bondi energy-momentum four-vector and show that indeed it is a 
Lorentzian four-vector. The details are given in appendix 7.2. 

In the discussion section, Sec. V, we summarize and discuss our results, 
including the results for the Einstein-Maxwell case. We also speculate on 
possible consequences of these results. One of our more interesting results - 
worth speculating on - is that a simple geometric condition on the Einstein- 
Maxwell fields leads to the Dirac value of the gyromagnetic ratio. 

Our method for the study of the physical identifications is based on the use 
of spherical harmonic expansions so that the coefficients in these expansions 
become identified with physical objects. It is clear from the non-linearity of the 
theory that it is impossible to work with the exact infinite expansions and, thus, 
we must work with a truncated series. Specifically we expand everything up to 
and including the I — 2 harmonics and then include only 2 nd order products. 
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(In one place, where the result is physically so attractive, we have included a 
3 rd order term.) Almost all of the non-linear terms in our results arise from the 
(frequent) use of Clebsch-Gordon expansions of spherical harmonic products. 



2 Preliminaries 
2.1 Null Infinity 

After integrating [51 OH ED [H] the Einstein (or Einstein-Maxwell) equations 
along null surfaces in the spin-coefficient formalism for large affine parameter, 
r, we are left, in the limit of large r with several evolution equations[T2] (for the 
asymptotic Weyl tensor components) that are referred to as the "asymptotic 
Bianchi identities" . (These equations are analogous to the remaining Maxwell 
equations[T2] after all the radial integrations have been performed.) This limit 
of large r has been formalized by Penrose into the idea of the future null bound- 
ary of space-time, a null three-surface, and referred to as 3 + . Our remaining 
equations are thus a set of differential equations on 3 + . This null boundary has 
the topology of S 2 x R and an assignment of (Bondi) coordinates (u, £, £ ), with 
u on the R part and (C, C) as complex stereographic coordinates on the S 2 . The 
freedom in the choice of these Bondi coordinates, known as the BMS group is 
geometrically the asymptotic symmetry group. The BMS group is composed of 
two parts, the supertranslations 



u = u + a((,() (1) 

(OC) = (00 

with a(£, £) an arbitrary smooth function on the sphere and the Lorentz transformations [H 
given by 
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Ku (2) 

i + OC 

« + b) (oC + b) + (cC + d) (cC + d) 

C = , , , i ad- be = l. 
cQ + d 

Though we are dealing with real space-times and 3 + is a real three-surface, 
since we assume that all the relevant functions are analytic, it is useful to allow 
the u to take on complex values close to the real and to allow £ to deviate 
slightly from the complex conjugate of £. 

In addition to the Bondi coordinates we also have a null tetrad (l a , n a , m a ,m a ) 
that is associated with the Bondi coordinates (and does change with a BMS co- 
ordinate transformation) where n is tangent to the null generators of 3 + , m a 
and m a are tangent to the u = constant slices' of 3 + while / is the null vector 
normal to the slices pointing into the space-time along the null surface u = 
constant. 
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2.2 Further Structures 



There are several important structures associated with J + that we now describe. 

• In addition to the Bondi 'slicing' of 3 + given by the 'cuts', u = constant, 
one can take an arbitrary real one-parameter, s, family of slices given by real 
cut functions, 

u = G(s,C,C), (3) 
or their generalization to analytic complex cut functions 

u = X(tXX) (4) 

with inverse functions 

r = T(u,CX) 

where the complex r must be able to be chosen so that u is real. The freedom 
in the choice of the parameter r: r = F(t) with F analytic, will be used later 
to normalize a physical variable. 

• At every point of 3 + there is the past light-cone of rays going back into 
the interior. The directions are labeled by the complex stereographic angle 
(L, L) with the zero value taken along I and the infinity along n. Any complex 
angle field L((, £), i.e., a stereographic angle given for each point on J + , can be 
expressed [UJ in terms of a complex cut function, restricted to real u's by 

L(u,C,C)=3(t)*(t,C,C) (5) 
t = T(uX,(). (6) 

The subscript (r) means the application of S holding r constant. 

• We will have considerable use for the local Lorentz transformation (null 
rotations) at each point of 3 + parametrized by the arbitrary angle field 

L = L(u,CX) 
that preserve the vector n, namely 



I* = I + Lrn + Lm + LLn (7) 
m* = m + Ln 
n* = n 

For each point on 3 + , the null vector I* determines a null geodesic extending 
backwards into the space-time so that the field of l*'s determines a null geodesic 
congruence of the space-time. 
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2.3 What lives on 3+? 

The following functions are defined [T^j on 3 + ; 

• it = o(u, £, C)i the asymptotic shear of the null geodesic congruence with 
Bondi tangent vector I. 

• tpi (u, (, () and -02 ( w i C) are the leading terms of two tetrad components 
of the Weyl tensor, ip 1 — 0° /r 4 +.... and tp 2 — V*§ A 3 + ■■• 

with 



V'i = -C abcd l a m b l c m d 
tp 2 = -C abcd M a n b l c m d 

Under the tetrad transformation, Eq.([7]), the Weyl components and 
transform as 

= V? - 2,Lipl + 3L 2 g(a)- + L 3 a" (8) 
= V^2 - 2LB{ay - L 2 W (9) 

The ip® and ?/>2 satisfy the asymptotic Bianchi Identities 

(^y =-^° 2 +2aB(ay (10) 

{i>ly -a{a)~ (11) 

^2 -^2 =9 2 <7-g 2 a + (cr)ff-(?)-(7. (12) 

The last two equations can be rewritten in terms of the 'mass aspect ', 

as 



*=o-ct, (13) 
= y = ijj% + d 2 7f + <r{ay. (14) 

When a Maxwell field is present, these equations become modified and the 
Maxwell equations must be included [12] : 

/,0v _ a„/,0 1 1 i;„,aOT 



)■ = -5^2 + 2ff3(ff)- + 2fc0 l 1 V2 (15) 

(V 2 °)- =-S 2 (a)-- ( x(a)-+fc^2 (16) 

* = * = 0° + B 2 a + a{wy, (17) 

fe = 2Gc~ 4 (18) 

+ 90° - 0$ = 0, (19) 

+ a$ = 0, (20) 
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2.4 Asymptotically Shear-free Null Geodesic Congruences 

• If the angle field L(u, £, () satisfies the differential equation[T51 [T7] 

BL + LL = a(u,(,() (21) 

the null geodesic congruence determined by the null vector field given by Eq. ([7]) 
is asymptotically shear-free. It has been shown earlier [T7] that solutions to 
Eq. (|2"Tj) that are regular on 3 + , i.e., have no infinities, are given by the following 
construction: 

L(u, C) is given parametrically by 

i(u,C,C) =9( t) X(t,C,C) (22) 
u = X(t,C,0 (23) 

where X(t, (, C) is found by first solving the 'good cut' equation [19]. 

d 2 Z = <j(Z,C,Q (24) 

whose solutions are known to depend on four arbitrary complex parameters, 
z a , i.e., Z = Z(z a ,(,C)- By choosing an arbitrary world-line in the parameter 
space, (ii-spacepj)]), i.e., z a = £ a (r), X(t,(,() is determined by 

u = X(T,(,0=Z(e(T),(,t)- (25) 

We thus have that every regular solution to the asymptotically shear-free con- 
dition, Eq. ( I2ip , is determined by an arbitrary complex world-line in a four- 
complex dimensional parameter space. The freedom, mentioned earlier, in the 
choice of t is used later to give a normalization to v a (r) = d T ^ a . 

We could reverse the procedure just described and assume that u = X(t,(,() 
was given rather than a(u, £)■ Then a(u, CiC) an d L(u, £) could be deter- 
mined parametrically by 

u = X(r,C,C) (26) 
L(u,C,C)=3 (t) X(t,(,0 (27) 

a( U ,C,C)=3(r)*(T,C,C) (28) 

or more explicitly by the harmonic series, 

u = X(rXX) = ±f{r) ~ \iKr)Y^,0 + e j {r)Y^,0 + 

(29) 

L(u,C,C) = C (r)^(CC) ~ 6C«(r)y a ^(C,C) + - (30) 
a(u,C,0 = 24e j (r)Y 2 %+..., (31) 

with £ a = (£°(r), £ l (r)).Tlie complex parameter r is given by r = T(u,(,0, 
Eq.(j6|). Note the important point that when r is replaced by T, the spherical 
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harmonic decomposition of those variables becomes non-trivial since it involves 
products of different spherical harmonics. 

In the special case of flat-space, (with cr{u, 00 = 0), the asymptotically 
shear-free congruences become shear-free congruences and the regular solutions 
to Eq. pTj) . given parametrically, become 



u = x(t,c,o = r(T)fa(c,o = -^e(r) - ^m^kco (32) 

L(u, C,C) = 3X(r, 00 = C(r)m a (C, C), (33) 
^ = ^(l,^,-^i,±±% (34) 

vs s; 2 v i + cc l + CC i + CC 

m Q (CO=3^(CO = ^(0,l-C 2 ,-*(l + C 2 ), 2C), (35) 

with l a (C, C) a flat-space null vector on the 'light-cone' composed of the I = 0, 1 
harmonics. In other words the regular flat-space shear- free null geodesic congru- 
ences are determined by an analytic complex curve in complex Minkowski space. 
One can reinterpret the Minkowski space curve as a curve in the space of asymp- 
totic complex Poincare translations. The advantage of this latter interpretation 
is that it applies just as well to the asymptotically flat space-times. 



2.5 Spherical Harmonic Decomposition 

All the functions on 3 + that we are dealing with have a spin weight (s) and 
most have a definite conformal weight (w). They thus can be expanded in the 
spin- weighted tensor harmonics- [20] . Denoting the spin and conformal weights 
by (s,w) as a subscript e.g., Wr S:W ), we have 



u = X ee X (0;1) = -^°(t) - if (r)yO(C + (r)l& (0 + - (36) 



a = a i2 ^ 2) =2ie j (r)Y 2 %+... (37) 

l ee l (1jX) = e'MnUoc) - 6f»y 2 yco + - (38) 

€ = €(o,-3) = T + + i;° 2ij Y 2 % + ... (39) 

< ^ ^(1,-3) = VW< + €i^2ii + - (40) 

* ee * (0 3) = *° + + *« y o + ... (41) 



^ = <-^+C y 2« + - (42) 
^ = Q + te° i +^ 2 y + --- (43) 

$ = &V+<i y a 7/ + »- (44) 
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The basic idea is to try to give physical meaning or significance to the har- 
monic coefficients. An important well known example of this is Bondi's identifi- 
cation of the I = (0, 1) parts of the mass aspect with the total energy-momentum 
four-vector of the interior sources, (Mc,P l ) which is explicitly given by 

* = _*/plM + (45) 



From the I = harmonic coefficient in Eg. ([13]) one obtains the Bondi mass 
loss equation which allows us to identify the coefficient as proportional 

to the 2 nd derivative of the mass quadrupole. 

A few other identifications coming from Maxwell theory are: Q is the 
Coulomb charge, while is proportional to the complex electromagnetic 

dipole moment (electric + i magnetic), 0^ is proportional to the 2 nd time 
derivative of ^ while <\>\^ is proportional to the 3 rd derivative of the quadrupole 
moment. 

Very roughly speaking we make the approximate identification of -0^, (which 
is complex) with the complex gravitational dipole moment, i.e., with the mass 
dipole plus l V times the angular momentum. 



3 Identifications 

Our starting point for the physical identifications is Ea.lfTO]) 

(V>?y = -3^ + 2<r3(^)- (46) 

and the spherical harmonic decomposition. We noted that the I = 1 part of ify 
was proportional to the Bondi three-momentum, P. The linearized version of 
Eq.(|3fi)) leads to the fact that the I = 1 part of (V'lj)' is also proportional to P. 
This in turn suggests that i/j^ itself should be, at least in the linearized version, 
proportional to the mass-dipole, MX, where X should be identified with the 
position vector of center of mass so that MX ' = P. 

This is our fundamental observation. It will be analyzed and generalized. 

We begin by postulating three different models or methods for the identi- 
fication. Though we have one method (the third) that we consider to be fun- 
damental and correct, nevertheless we felt that at least a few others should be 
explored to see if they would be reasonable choices. Our criteria for selecting a 
model are: (i) it should predict already known laws or reasonable new laws, (ii) 
it should have a clear geometric foundation and have a logical consistency and 
(iii) it should agree with special cases, mainly the algebraically special metrics 
or analogies with flat-space Maxwell theory 

The first (naive) approach is to simply assume that the Bondi 3-momcntum 
is given by P — MX ' = M V . The second approach is based on the flat-space 
transformation (translation) properties of dipoles, i.e., transforming to the cen- 
ter of mass or charge, while the third approach generalizes this to asymptotically 
flat space-times. The first two do not satisfy our three criteria while the third 
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one does. We nevertheless felt it was worthwhile to see how the increasing 
sophistication of the methods led to improved results. 

Remark 1. We have (for notational simplicity) totally abused standard 
notation. We allow the indices (i,j,k,l...), which are Euclidean, to be raised 
and lowered with impunity and allow repeated indices, e.g., v k t; tk , to indicate 
summation. 

Remark 2. Though u is the conventional Bondi time, it is more appropriate 
to use w — \/2u, it being the retarded time. Derivatives with respect to u are 
denoted by dot, ('), while w derivatives are given by a prime, ('), i.e., (')= \/2(')- 

Claim A very important computational issue is to find the inversion of 
Eg. \3b\ ). i.e., to find t = T(u, CO- The reason for its importance is that we 
must be able to explicitly eliminate the r and replace it by the Bondi u (or w) in 
the parametric expressions for L(u, £), i- e -j in the expression for the complex 
world-line £ a (r). The approximate inversion (linear) is given J21^by 

t = T(u,(,0 = w + H^(CC) - V2^{w)Y^(U) + (47) 
3.1 The First Identification Method 

For our first (toy) model, we assume that the Bondi momentum (at least at 
linear order) has the standard kinematic form, namely 

P = MV, (48) 

and try to find the consequences. From our perspective this already is unsat- 
isfactory, since we would like this kinematic expression to follow directly from 
the Einstein equations rather than to assume it. In fact, in the second and 
third models this relationship is a derived result. Nevertheless from a heuris- 
tic point of view we believe it worthwhile to show the line of reasoning — in 
this simpler case, - before is used later for the preferred models. In addition, 
roughly speaking, from this first model we can see what consequences can be 
anticipated. 

Thus from the above argument, we take ipu to have the form 

ip% = aMXi (49) 

where a is a constant to be determined, M is the Bondi mass and Xi is a complex 
three- vector in some unknown space (hopefully) to be determined. The real part 
of Xi is tentatively associated with the 'center of mass position' and its derivative 
with the 'center of mass' velocity, i.e., X' Ri — Vi =>- V . 

Using the harmonic expansions, Eqs. ([39|) . (|40|) and Clebsch-Gordon prod- 
ucts, [see appendix C], Eq. (|46| becomes 

A = ^-^r'V m e yfc . (50) 
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while the expansion of the mass aspect, Eq. ljlTj) , is 

* = TH + ^fH^f' + + ^^^f'eu^l + ... (51) 

Using Eqs.(??), ([50|) and (|49|) we identify, from Eq. (j45|) . the mass and three- 
momentum 



2V2G 5 G 5 

_k , r , / .192 i—jrnf.jr,-, 

P k = MX' k + i—c 3 e £ lm e ljk 



a 
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Writing the variables as 



T = T R + iT 

^k = ^Rk + i^ik 
r-lm elm i ■ Am 

? = €r +Kr 

gml = v lm = v lni _ 



the reality conditions, Eq. (fT4]) . becomes 

M = M -^^($d? (52) 

Ti= (2^2 { ^_^ vf) (53) 

MX' Ik = -Q^if&vt + &v?)e ljk (55) 



Finally, Eq-dJ, 



(^y = -g 2 (cf)--^r, 



after using the Clebsch-Gordon expansion for cr(a) ", yields for the I = 0, 1 terms, 
the mass and momentum loss expressions, 

M '=-lf(«+«) (56) 
pfc/ = pk _ l^ (w « w « _ „« t ,«)e (Ifc . (57) 
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e iJ = (a + itf ) = T^rAQZL + iQXiJ- (59) 



There are now two identifications that can be made immediately from these 
expressions. From the quadrupolc energy loss relations[25j 

Q 

M = —^-y{QMassQMass + Q SpmQ Spin) (58) 

we can relate to the mass and spin quadrupole moments by 

G_ 

This identification applies to all three models. 

Furthermore, we can substitute the expression for P. i.e., Eq. (|54[) , into 
the momentum loss equation, Eq. (|57[) . obtaining equations of motion for the 
'position' vector A^, 

(04^2 2 1Q9r 2 

MX'ku ~ K -^(i? l vT - &vr)'e lik = ^-(v l /v% - v l lv?)e ilk . (60) 

This can be interpreted as Newton's 2nd law, with the second term on the 
left coming from gravitational radiation reaction and the right-side from the 
'rocket' or recoil force due to momentum loss. 

For the last identification in this model we consider Eq. |55|) . From the 
analogy with the rotating mass solution, (with and without charge [7]), and 
results, which will be described later, we identify 

S k = McXik 

as the spin angular momentum so that Eq. l|55[) becomes the angular momentum 
loss equation, 

Summarizing the results for this physical model; we have assumed that the 
(complex) I = 1 part of ip® k was proportional to the mass dipole moment and 
the spin angular momentum. The Bianchi identity Eq.( l46p allowed us to find a 
kinematic expression for the Bondi linear momentum in terms of the mass dipole 
moment (or center of mass) while the Bianchi identity, Eq.f ITT]) . yielded the 
equations of motion. The I — 1 reality condition became the angular momentum 
loss equation. 

Though this model does lead to physical identifications that are quite rea- 
sonable, it nevertheless contains several very unsatisfactory features. First of 
all there is no hint of a geometric meaning for the 'position' vector Aj,. What 
space is it in? It was pulled out of thin air, just so that the momentum could 
be written as P = M V . The constant a was arbitrary until determined by the 
kinematic meaning of P. A second deficiency is that the angular momentum 
is given only by the spin without any suggestion of an orbital part. Lastly, 
this identification is at total odds with the identifications obtained for the alge- 
braically special metrics 26, 28, 2TJ [2TJ or with the analogous flat space Maxwell 
identification with the center of charge[6]. 

We will sec that the third model docs correct these deficiencies. 
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3.2 The Second Identification Method 



Our second and third identification models are based on complex curves, in the 
second model these curves are associated with flat-space shear-free null geodesic 
congruences while in the third model they are associated with the asymptotically 
shear-free null geodesic congruences of the relevant space-time. These curves, 
as we mentioned earlier, are geometrically well defined as complex world-lines 
in the space of complex Poincare translations at 3 + . The ideas involved are 
not obvious and do lie outside of the usual default understandings and thus 
it takes a bit of 'getting used to'. Basically, the idea is to generalize the 
standard action of moving the origin of coordinates to a new point, (the center 
of mass), so that the mass dipole vanishes for that origin. The generalization 
consists of making these transformations (translations) complex and moving to 
the complex center of mass, so that not only does the mass dipole vanish there, 
but so does the angular momentum. This type of procedure, which worked 
perfectly, in linearized gravity and Maxwell theory [22, 23J must now be applied 
to asymptotically flat space-times. 

We begin with the observation that most attempts [2] at defining the angular 
momentum have the imaginary part of the I = 1 harmonic of V'? as part of the 
definition. 

Remark 3. Most of the angular momentum definitions based on symmetry 
arguments start with the 1=1 part of ip® and augment it with terms quadratic 
in the shear and its derivatives. Later in this section we will return to the issue 
of these quadratic terms. See Appendix 7.4- 

The central part of our argument is to use the I = 1 part of the null rotated 
Weyl component, ; from Eq.©, i.e., 

= V? - 3L^° + 3L 2 V° - L 3 ip1 (61) 

as the transformed (translated) complex dipole moment and then set the I = 1 
part of -01* to zero to (partially) determine the L . (In flat space this type of 
operation really moves the origin of Bondi light-cones to a new complex origin in 
complex Minkowski space so that the center of charge and the magnetic dipole 
both vanish.) In other words, we will set tp°* = so that the I = 1 part of 
ipi is given by 

0° |, =1 = (3L0 2 ° - 3LV§ + L 3 i/>°)\ 1=1 = 3(L$)li=i~ (62) 

It is this relationship that had, in the first model, its counterpart in the arbitrary 
assumption of Eq.([49|. Now it has a justification and a geometric meaning. 
Also it does not contain an arbitrary factor a. The quantity ip® is uniquely 
defined by the choice of L(u, £, (). 

Our second and third models depend on two different, but related, choices of 
L(u, C, C)- For the third model, we choose the L{u, £, £) to determine an arbitrary 
asymptotically shear- free null geodesic congruence, via Eqs. (|22j) and (|25|) . while 
for the second model we choose the L(u, CiC) associated with flat-space shear- 
free null geodesic congruences, Eqs. (|33p and (|32p . The L{u, Q,Q in either case 
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is determined by an arbitrary complex world- line £ a (r), which contains the 
information about the center of mass world-line and the spin. Though the third 
model is more consistent and has a more logical basis, for two reasons we first 
discuss model two. It is easier to work with and the final results are very similar 
to those of the third model. 

Using the flat-space L(u, £, £), working with terms up to 2 nd order, with 
frequent use of Clebsch- Gordon expansions, we find, from Eq. (|62"]) , that 

• 1 ^ , .3V2 , n A 18 

5 



V.-?, = 3T[f H + i-e kji vH j ] + ~ ■ (63) 



There are two things to immediately notice: 1. since T is proportional to the 
mass, the first term has exactly the same form as in the first model with £* now 
replacing the X 1 as the 'complex center of mass' and 2. the ip^ is fully determined 
without the constant a. In addition, new terms have appeared, some of which 
will have immediate physical meaning. From the Bianchi Identity, Eq. (|10p . 
using Eq. (|63p and the reality conditions, we have for the I = component of 

M = Mo-^(^+^f) 

_ 24(24)V2 , jj jj, {-ij+ijis 
oc 

which are identical to those of the first model. However for the I = 1 components 
new terms appear: 

For the momentum P we have 

P k = M v% + - - (e R r4)%jk + R\ (64) 

R k = + $&Y - ~ (65) 

The / = 1 part of reality condition yields an expression that consists of two types 
of terms, total time derivatives and the others. We define it as the conservation 
of angular momentum equation: 

j kl = + $&%ik, (66) 

J k = M cff + M (44 " CA)e ijk - - m k )■ (67) 

The angular momentum flux is the same as in the first model. However the 
angular momentum itself has much more physical content. The first term 



15 



(68) 



is again identified as the intrinsic spin, while the second term is precisely the 
orbital angular momentum, r x p, while the third term is the contribution 
from spin-precession. The fourth term is new, involving dipole and quadrupole 
coupling. This could be considered as a prediction of the theory though it 
probably is untestable. With these identifications, going back to P k in Eq. (j6"4"|) , 
we observe that the second term 

^'^ijk = r 2 J = C- 2 (S'.w) fc 

is the Mathisson-Papapetrou contribution to the linear momentum. 

Remark 4. Many of the kinematic expressions for the physical quantities 
have been given up to second order. When the dynamic equations are used, it 
often turns out that these quantities are really of higher order and should be 
neglected for second order considerations. We felt that for purposes of under- 
standing their kinematic role, their inclusion here was important. 

The equations of motion for the position vector, are obtained by substi- 
tuting the P from Eq. ([64j) into the momentum loss equation 

Phi T^k 192c , ij ij [j ij . 

= F = -^r( V I V R ~ v R v I > e ilk 

yielding 

Mov% ~ ^(&vf + + R k ' = ^vMatk- 

This again is Newton's 2 nd law but with spin-coupling forces, radiation reaction 
and momentum recoil. 

Though this second model is far superior to the first model it still has certain 
deficiencies. It does not yield the correct dynamics for the algebraically special 
metrics. In other words the dynamical equations would differ from those of 
the algebraically special metrics. Also there is a lack of logical consistency in 
this treatment. We have used the L(u, £, £) for the flat-space shear-free null 
geodesic congruence. We should have used the L(u, (, £) associated with the 
asymptotically shear-free congruence, i.e., the L(u, £, £) from Eq.(l22jl. We would 
then have the null rotated Weyl component, ipii? based on an asymptotically 
shear-free null congruence so that the quadratic terms in the shear, that would 
normally appear in the definition of angular momentum, would be now absent, 
[See Remark 3.] and our procedure of setting i/j^ = would be consistent with 
all definitions of center of mass and angular momentum at the complex 'origin'. 

3.3 The Third Identification Method or Physical Model 

The third method is, as we mentioned earlier, basically the same as the second 
method. So rather than repeat ourselves with the vacuum ' asymptotically shear- 
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free', L(u,C,0> we w ^ S^ ve the results for the asymptotically-flat Einstein- 
Maxwell equations. Since the method and arguments are so similar to those of 
method two and the detailed calculations are so long, we will simply summarize 
the results. 

The physical identifications are substantially changed by the presence of the 
Maxwell field. For the vacuum case our results apply to the general asymptot- 
ically flat situation. For asymptotically flat GR, with the Maxwell field, the 
situation is different, ft turns out that the Weyl tensor condition (f6"2")) has a 
counterpart for the Maxwell tensor. In other words, the / = 1 part of a null 
rotated <Pq, i.e., </>§* when set to zero, determines a different complex world-line 
(the complex center of charge) for which the electric and magnetic dipole mo- 
ments vanish. We consider only the special case where the gravitational and 
electromagnetic world-lines coincide. In some sense this coincidence of world- 
lines implies that the source has a restricted structure and is relatively 'simple'. 
We take this as the meaning of a gravitational elementary particle. 

We start with a very general cut function and its inverse that came from a 
given Bondi shear, i.e., from Eqs.([23) and 

w „, - s 1 j,, , 1 



u = ^= = X(t, C, C) = ^°(r) ~ ^to 1 "^ + ^» y 2°y(C, 



(69) 



r = T(u, C, C) = w + ^fH^K, C) - V2^(w)Y 1 %((, C) + ... (70) 



V2 
2 

l = ew^co - 6fw 2 yc,o (71) 

C = (f, €(t) = (r,f (r) = Cr(w) + it?(w) (72) 
C' = v a (w)=v%{w)+iv<}{w) (73) 

After eliminating r, Eq. ([62]) yields the expression for tjj^ : 



= 3T[C H + i\e kji v k e + + i^ eiji ip° 2j e - H^o..^' (74) 



5 %fcj? 5 V2jK 



Ni = 6^^ W _ W2 /; ^ + i™l e . miV V C k (75) 

Aside: Later, for the vacuum case, we will associate the imaginary part 
of Eq. ([7^|) as the total angular momentum. 

Integrating the Maxwell equations leads to 
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i 2 72 

<t>° u = V2Q[v\w) + i-eqiv*? + N if ] + iv^Qey,^* - - yQwV' 

(76) 

, . 24-^2 „ , fc ,24\/2 , ■ , 
+ l—^r-ZljWokjZ ~ miVOmlS ) 

00. = _2Q[„* + i-eyjW)' + ^V 1 "] - £2Qe fcii (^ fe )' + ^(^ fc )" 

+ ^Q(tA,« )' - gie w (^C tt )' + if e, K ($U im )" 
From the Bianchi Identity Eg. ([15]) we obtain 



P k = mv\ + ^(v^ - - (^yy m - |f < (77) 



c 



+ %\ 2 &r - M' + v^A'^k + n fc (78) 



15V2c * Si SK/ lJK 5c v * ^ ' J iU ; ' 5c 1 

4 (rs ii "T) il " 1 n iJ "n u '"v,-, 36V2c 2 , 

45c 4^M -E i e bfe 5( -f 1 
2(24) 2 C 2 , ij c u 



and -D^ are respectively the electric quadrupole and magnetic quadrupole 
moments found from the I — 2 radiation term in the solution of the Maxwell 
equations. All non-linear terms involving the quadrupole terms are gathered 
into the H k . 

The vanishing of the imaginary part of the reality condition yields the rela- 
tions 

(79) 
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where J k , identified (from the dynamics rather than through the conventional 
symmetry argument) as the total angular momentum, is given by 

J k = Mc£j + M(&4 - Q4)e ijk + ^v k - (ft< + ^}4')e ijk + K k 

(80) 



Y^i (ft<"' + e>^'" + " + *V2(2&D%'" - 2#r>«'" 



36V2 



.3 



The justification for calling this the angular momentum is the same as in the 
previous model, coming from the sum of terms interpretable as the spin, orbital 
and precessional moments. 

In the absence of a Maxwell field J k is simply proportional to the imaginary 
part of 

V2c 3 
12G 

Remark 5. In the published literature^ there are ambiguities in the def- 
inition of angular momentum. For us, because of our approximations, these 
ambiguities disappear, i.e., the ambiguous terms are higher order. See Appendix 
74- ' 

Our last results arise from are the mass/energy loss equation and the mo- 
mentum loss equation. The later being the dynamic equations for £? R . By sub- 
stituting the P k of Eq. (|77|) into the Bianchi identity, Eq. (fT6|) we find 



J l = -^77^u\l- (81) 



= - ^{QkassQMass + Q SpinQ Spin) ~ ^^ V R V R + V I V I ) ( 82 ) 



■ \ e e m i 



180c 7 

F k = 2G (QljHI Qijlll _ QljUl Qijlll N _ Q_, U A, _ Jl,,H\ f . u 

_ 15c 6 Spi " aSS ^Mass^Spin) e ^k 3 4 V I R V R V I) e ilk 

(83) 



-(v j, D jk '" 4- v J/ D jk "') + —^—!D lj "'D ijl " - D lj, "D ij '")e n 



15c 5 v " ^ J ; ' 540c 6 

where we have now replaced the £ y by the more physical variables, the quadrupoles, 
>-/" 

G 



(Qtlass > Qlspin ) via 
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Writing out the equations of motion in detail is long and not completely 
enlightening, There are many non-linear terms whose meanings, other than 
they are interpretable as gravitational radiation reaction, are not clear. Instead 
we will write out a truncated version of Eq. ([55]) , hiding these terms in a single 
symbol, R k , so that we have 

Mv% + v\M' - ^-4" + R k = F k . (85) 

There are several things to notice here. First of all strictly speaking, we should 
ignore the term v R M' since it is third order because M' already is second 
order from Eq. (|82|) . We keep it, understanding its suspect nature, because it is 
potentially so important. For the moment ignore i? fe and F k and consider only 
the second term in the mass loss, i.e., the electric dipole energy loss, so that 
Eq.([8"5")) becomes 



Mv R - v R —rv R v R - —tv r = 0. (86) 



These[231 123 are the classical equations of motion for a charged particle that 
contain the well known electromagnetic radiation reaction force and exhibit the 
unstable run-away behavior. (The cubic term, though not sufficient large to sta- 
bilize the equation, has the correct sign for stabilization.) This classical equation 
has been obtained without the usual model building or mass renormalization. 
Returning to Eq. (|85|) , we see that in addition to the i? fc and F k , there are the 
extra terms in M' coming from gravitational and electromagnetic quadrupole 
radiation. It is very hard to see the consequences these terms though one can 
see that they are on the side of stabilization. It is then easy to conjecture 
that coupling electrodynamics with general relativity stabilizes the equations of 
motion. 

There are other comments concerning our results that should be made. 

• The mass/energy loss equation contains both the gravitational radiation 
quadrupole expression (this has been adjusted by definition) and the classical 
electromagnetic dipole and quadrupole radiation expressions that come straight 
from the construction. This allows us to identify 

(i k = Qi k (87) 

as the magnetic dipole moment, so that with the spin definition, S k — Mc^j, 
we obtain the Dirac value of the gyromagnetic ratio, i.e., 

g = 2. (88) 

• The angular momentum expression, Eq. ()80|) . contains the spin, the orbital 
angular momentum and a precession contribution. There is now a prediction 
that there is a charge-spin contribution to the total J, i.e., j^jS k ' , as well has 
higher order corrections. 

• Probably the strongest argument for the validity of this approach to ex- 
tracting physical information from the asymptotic field is the observation that 
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in the angular momentum loss equation, Eq. ([75]) , the dipole contribution to the 
angular momentum flux coming just from the electromagnetic field, i.e., 



^L.( v ^ R+ v}'4)e ljk (89) 

exactly coincides with the classical electrodynamic angular momentum flux |25j . 
This result was obtained with no a priori expectations. 



4 BMS Invariance 

Our results concerning the identifications of physical quantities from the asymp- 
totic fields were all obtained in an arbitrary but specific Bondi coordinate/tetrad 
system. The question is what are the relations between the same quantities but 
calculated in a different Bondi coordinate/tetrad system? In other words, we 
want to know the transformation properties of our physical variables under the 
action of the BMS group[l EJ [15l [TBI El]- As we pointed out in Sec. II, the 
BMS group is composed of two parts, the supertranslations and the Lorentz 
transformations given respectively by 

u = u + a((,C) (90) 

(OC) = (00 

with a(£,0 an arbitrary smooth function on the sphere considered now to be 
small i.e., as a first order quantity and 

u = Ku (91) 

K= l + ^ „ 

« + b) (oC + b) + (cC + d) (cC + d) 

C = —r— ad-bc=l. 
cC + a 

with (a, 6, c, d) the complex parameters of SL(2,C). The invariance under super- 
translations is actually subtle and requires a bit of thought while the invariance 
under the Lorentz group, though straightforward, requires some more techni- 
cal background. Since all our calculations were done under the assumption of 
second order perturbations off a Reissner-Nordstrom background, we must keep 
the BMS transformations small, i.e., close to the identity. 

4.1 The Supertranslations 

We first treat supertranslation invariance. 

Starting on 3 + with some arbitrary but given Bondi coordinates (u, £, C) an d 
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tetrad (I, n, m, m) we saw that there was a null rotation, Eq.0, 

1* = I + Lfn + ~Lm + LLn (92) 
m* = to + Ln 
n* = 77 

that was determined by the choice of the angle field L(u,£, £). The L(u, £,C) 
was then determined first by the asymptotic shear-free condition, leaving the 
freedom in the I = (0, 1) harmonic coefficients, and then by the requirement 
that the I = 1 harmonic component of the 'rotated' Weyl component, ip®* , 
i.e., ipi* = 0, should vanish. The important point to note is that the Weyl 
tensor component ip®* is determined by the new tetrad, (I* , n* , m* ,rn*) , It is a 
geometric structure given independent of the choice of coordinates - depending 
on the fixed known (*)-tetrad. If we have a second Bondi system with Bondi 
tetrad (I, n, m, m), it will have been obtained by a different null rotation via 
some other angle field, L, from the (I, n, m,m) tetrad. This simply means that 
there is a different null rotation, going now from the Q-tetrad to the previously 
determined (*)-tetrad. The ipi* |i=i thus remains zero when this harmonic is 
extracted from the ip®* when the 'u' is held constant. There however is a serious 
issue that must be raised. When we extracted the I = 1 part of , it was done 
at constant value of the Bondi l u' . But now, after the supertranslation, the 
harmonic decomposition should be done at constant t u'. However, it turns out, 
because of the first order BMS supertranslation, u = u + a(C,C)> the ip°* |;=i 
part of ipi* , now obtained by holding 'u' constant, remains zero up to second 
order. So up to our accuracy, our results are supertranslation invariant. The 
more difficult issue of how to deal with finite supertranslations will be discussed 
in the conclusion. 

4.2 Lorentz Transformations 

To deal with the Lorentz subgroup of the BMS group requires a review of the 
representation theory. 

The theory of the representations of the Lorentz group was beautifully de- 
scribed by Gelfand, Graev and Vilenkin|13j using homogeneous functions of two 
complex variables as the representation space. We will summarize these ideas 
using an equivalent method 14 [ namely by using spin-weighted functions of the 
sphere as the representation spaces. Representations are labeled by either of two 
numbers (ni, 112) or (s, w), with (771, 77,2) = (w—s+1, w+s+l). The 's' is referred 
to as the spin weight and 'w J as the conformal weight. The representations are 
referred to as D( ni n2 y The special case of irreducible unitary representations, 
which occur when (ni, 772) are not integers, are not of interest to us and will not 
be discussed. We will consider only the case when (771,772) are integers so that 
the (s, w) then can take on integer or half integer values. The representation 
space, for each (s,w), are the functions on the sphere, ?7( Sltu )(C) C)i that can be 
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expanded in spin-weighted spherical harmonics, s Yim(C, 0) so that 

oo 

sYUCX), (93) 
and transform under the Lorentz group, Eq. (|9ip . as 

V(s, w )(lC)^e lsX K w V ^ w) (CX) (94) 

AX c C + d 



with 



K = (l + Cm< + b)(aC + b) + (cC + d)(cC + d)}- 1 , 



(95) 



These representations, in general, are neither irreducible nor totally reducible. 
For us the important point is that many of these representations do possess an in- 
variant finite-dimensional subspace which (often) corresponds to the usual finite 
dimensional tensor representation space. Under the transformation. (|M|) . the 
finite number of coefficients in these subspaces transform among themselves. It 
is this fact which we heavily utilize. More specifically we have two related situa- 
tions: (1) when the (n 1; n 2 ) are both positive integers, (or w > |s|), there will be 
a finite dimensional invariant subspace, Df and (2) when the (— m, —n 2 ) 
are both negative integers there will be an infinite dimensional invariant sub- 
space, D7_ ni _ n2 v One, however, can obtain a finite dimensional representation 
for each negative integer case by the following construction: One forms the fac- 
tor space, Z?(_ nii _ n2 -j/Z)7 ni _„ 2 y This space is isomorphic to one of the finite 
dimensional spaces associated with the positive integers. The explicit form of 
the isomorphism, which is not needed here, is given in Held & Newman 14 . 

Rather than give the full description of these invariant subspaces, which is 
available elsewhere, we will confine ourselves to the few cases of relevance to us. 

I. For s — and w — 1, [(ni,n2) = (2,2)], the harmonics, / = (0,1) form 
the invariant subspace. Applying this to the cut function, X(t, £,£),we obtain 
from 



« = *<o,u = ^V) - \e(r)Y°(<:,Z) +e 3 (T)Y 2 %((,o + ... 

(96) 

Invar, subspace = ^V) - (r)y£(£C) (97) 

(r) = (£° (t) , f (r) ) = Lorentz vector (98) 

This allows us to single out, in a Lorentz invariant manner, the four I — (0, 1) 
harmonic coefficients of the cut-function X(t, £, £) as a complex position vector. 
II. The mass aspect, 

* = * (0f _ 3) = *° + VYft + >I'"Vj;, + ... (99) 
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is an s = and w — —3, [(711,712) = (—2,-2)] quantity. The factor space is 
isomorphic to the finite dimensional positive integer space, [(711,712) = (2,2)] 
and hence the harmonics coefficients of I = (0, 1) lie in the invariant subspace. 
From the isomorphism, (which does change the numerical coefficients) we can 
construct functions of the form Eq. (|97l ) , which, in turn, lead to 

P a = (Mc, P*) = Lorcntz vector 

This gives the justification for calling the I = (0, 1) harmonics of the mass 
aspect, a Lorentzian four- vector, P a . 

III. The Weyl tensor component, has s = 1 and w — —3, [(7ii,7i 2 ) = 
(—3,-1)]. The associated finite dimensional factor space is isomorphic to the 
finite part of the s = — 1, w = 1, [(ii, 712) = (3, 1)] representation. We have that 

^1 = < ( i,-3) = + vujy-h + - ( 10 °) 

leads to the invariant subspace 

Invariant subspace = ipiiY-^ 1 . 

The question of what finite tensor transformation does this correspond to is 
slightly more complicated than that of the previous examples of Lorentzian vec- 
tors. In fact it corresponds to the Lorentz transformations applied to (complex) 
self-dual antisymmetric two-index tensors. As an example from Maxwell theory, 
from a given E and B, the Maxwell tensor, F ab , and then its self-dual version 

jyab+ _ pab _|_ j^p*ab 

can be constructed. A Lorentz transformation applied to W ah+ is equivalent |25j. 
(see appendix B) to the same transformation applied to 

V^t = (E + (101) 

These observations allow us to assign invariant physical meaning to our 
identifications of the position vector, £ a , the Bondi momentum, P a and the 
angular momentum J 1 . 

5 Discussion 

Starting from a very unorthodox point of view, we have tried to describe in 
the context of GR, (either in the vacuum case or for the Einstein-Maxwell equa- 
tions), the equations of motion of an isolated charged, massive body (our gravita- 
tional elementary particle) that possesses both an intrinsic spin and quadrupole 
moments and can radiate both gravitational and electromagnetic radiation. The 
point of view arises from considering a generalization of the algebraically special 
metrics where one can identify physical quantities and determine their evolution 
from the asymptotic field equations. 
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For the algebraically special metrics, the shear-free null geodesic congru- 
ence that are associated with these metrics, automatically assigned kinematic 
variables (e.g., a position vector and an intrinsic spin) to the Bondi energy- 
momentum four- vector. The Bondi evolution equations then became the equa- 
tions of motion. Our generalization consisted in observing that the existence 
of a shear-free null geodesic congruence could be generalized to the existence 
of an asymptotically shear-free null geodesic congruence. Applying the same 
physical identifications as with the shear-free congruences leads in exactly the 
same manner to the more general equations of motion. 

To test out if this result was accidental and, perhaps, other physical as- 
signments would lead to similar equations of motion, we tried two alternate 
strategies. They both were found to be lacking. They were not as natural or as 
physically meaningful as the third identification method. 

Since we are far from any of the standard or default approaches to the 
description of motion in GR it would be appropriate to summarize our results. 

Looking at the asymptotic Einstein-Maxwell equations, the asymptotic shear- 
free conditions lead, in general, to two different complex world-lines in the space 
of complex Poincare translations acting on 3 + , one from the Weyl tensor, the 
other from the Maxwell tensor. We have considered only the special case where 
the two world-lines coincide - defining this case as "elementary particles in GR" . 

Some results among others are: 

• The mass has a kinematic correction term dependent on the variable 
quadrupole moment. This could perhaps be considered as a prediction. 

•The Bondi linear three-momentum is expressed in kinematic variables, e.g., 
Mv k , ^s-Vji and the Mathisson-Papapetrou spin coupling, among others. 

• The imaginary part of the complex position vector is identified with the 
specific intrinsic spin angular momentum. From the solutions to Maxwell's 
equations the magnetic moment is seen to be the charge, q, times the imaginary 
part of the position vector. This agrees with the algebraically special charged 
spinning metric and leads to the Dirac value of the gyromagnetic ratio, i.e., 
g = 2. Though earlier we have defined a 'gravitational elementary particle' from 
this result, it should be noted that in the elementary particle community[29 it 
has been speculated that all charged elementary particles with spin have this 
property. One then has the question, what if any, is the relationship between 
the two types of particle? 

• One of the strongest arguments for our interpretations comes from the 
I = 1 part of the reality conditions, which is interpreted as the dynamics (the 
conservation law) for the total angular momentum. There is a total time deriva- 
tive term of a quantity, J k , that we define as the total angular momentum. It 
contains the spin, the orbital angular momentum and a precession term. The 
angular momentum flux contains three terms which come, respectively, from 
the gravitational quadrupole radiation, the electromagnetic quadrupole radia- 
tion and a term arising from the electromagnetic (electric and magnetic) dipole 
radiation. This latter term is identical to that calculated purely from electro- 
magnetic theory [25], 
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• From the Bondi mass loss equation, we can identify from the flux terms, 
the gravitational quadrupole but also see that our identification of the electro- 
magnetic dipole moments agrees with the predicted dipole energy loss. 

• From the Bondi momentum loss we obtain the equations of motion. In 
a sense we "derive" Newton's 2 nd law , F = Mv, where the force is a com- 
bination of electromagnetic radiation reaction [2H H5J, gravitational radiation 
reaction and a "rocket" recoil force from the electromagnetic and gravitational 
momentum loss. 

• Finally, we showed that each of the quantities that were identified as 
physical variables transformed appropriately under the Lorentz group, i.e., as 
Lorentzian tensorial objects. 

• Our results follow from the existence of a well-defined geometric structure, 
namely the UCF, a unique one-complex parameter family of slices of null infinity. 
This suggests that the higher order coefficients in the harmonic expansion of the 
UCF, e.g., £ lJ , ^ k , should be identified with time-derivatives of the higher 
multipole moments. 

There are other unfamiliar terms that could be thought of as predictions of 
this theoretical construct. How to possibly measure them is not at all clear. 

• One interesting physical prediction concerns the contribution that the 
charge makes to the total angular momentum. Looking at the equation defining 
J k , Eq.fJSO]), we see the (linear) contribution from 

which is a coupling between the charge and the changing magnetic dipole mo- 
ment. 

The final item to be discussed concerns the invariance under the supertrans- 
lation subgroup of the BMS group. 

An obvious question concerning the material described here concerns the 
issue of the extraction of the I = (0, 1) harmonics from different Weyl tensor 
components or from the universal cut function, X. We have consistently per- 
formed the extraction of the harmonic components on the cuts, u = const., or 
on neighboring cuts, u — const. + A, that are close to the first set of cuts. To 
2 nd order the results are unchanged. If we do go to arbitrary cuts, there is no 
reason for the extraction to lead to the same results. Our results thus appear 
to depend on the choice of cuts. In fact there is a canonical choice of cuts, i.e., 
a special one-parameter family of cuts, labeled by 's', on which the extraction 
should always be performed and for which there is no ambiguity. In the text, 
the u — const, cuts were sufficiently close to the canonical choice, so that, to 
2 nd order, they were the same. In fact, in principle, we should have been doing 
all our calculations on the s = const, cuts. It was, however, easier doing it with 
u = const. 

The question then is what are these canonical cuts? Returning to the com- 
plex universal cut-function 

u = X(t,C,0 
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we saw earlier that the r had to be chosen so that the u had real values. If we 
write 

T = s + iA(s,C,0 (102) 
then one can show {see appendix A} that A(s, (, Q can be chosen so that 

u = X( s + t A( s ,(X)XX) =X(sX,0 (103) 

is a real function of the real variable 's'. This is the construction of our canonical 
slicing. 

6 Acknowledgments 

G.S.O. acknowledges the financial support from CONACYT and Sistema Na- 
cional de Investigadores (SNI-Mexico) . C.K. thanks CONICET and SECY- 
TUNC for support. 

References 

[1] H. Wcyl, Space-Time-Matter, Methucn, London , (1922). 

[2] L. B. Szabados, "Quasi-Local Energy-Momentum and Angular Momentum 
in GR: A Review Article" , in Living Reviews in Relativity, June 2005 

[3] Adam D. Heifer, Gen. Relativity & Gravity, 39, 2125-2147, (2007). 

[4] Dray, T., and Streubel, M., Class. Quantum Grav., 1, 15-26, (1984). 

[5] Ashtekar, A., and Winicour, J., J. Math. Phys., 23, 2410-2417, (1982). 

[6] C. Kozameh, E.T. Newman and Gilberto Silva-Ortigoza, 
Class.Quantum.Grav., 24, 5479-5493, (2007). 

[7] E.T. Newman, E. Couch, K. Chinnapared, A. Exton, A. Prakash, R. Tor- 
rence, J.Math.Phys.6, 918, (1965). 

[8] H. Bondi, M.G.J, van der Burg and A.W.K. Metzner, Proc. R. Soc, London 
Scr. A 269, (1962) 21. 

[9] R. Sachs, Proc. Roy. Soc. 270, 103, (1962). 

[10] Ezra T. Newman and Roger Penrose, J.Math.Phys, 3 566-768, (1962). 

[11] E. T. Newman and T. Unti, J. Math. Phys. 3, 891, (1962). 

[12] E.T. Newman, K.P.Tod, in General Relativity and Gravitation, Vol 2, Alan 
Held, (ed.), Plenum (1980). 



27 



[13] I.M.Gelfand, M.I.Graev and N.Ya.Vilcnkin, in Generalized Functions, vol. 
5: Integral geometry and problems of representation theory, Academic Press 
, New York, (1966). 

[14] A. Held, E. T. Newman and R. Posadas, J. Math. Phys. 11, 3145 (1970). 
[15] R. Sachs, Proc. Roy. Soc. 270, 103, (1962). 
[16] R. Penrose, Phys. Rev. Ltt. 10, 66, (1963). 

[17] C. Kozamch and E.T. Newman, Class. Quantum Grav. 22, 4659-4665, 

(2005) . 

[18] B. Aronson and E. T. Newman, J. Math. Phys. 13, 1847-51, (1972). 

[19] R. Hansen, E.T. Newman, R. Penrose and K.P. Tod, Proc. Roy. Soc. A 
363, 445, (1978). 

[20] E.T. Newman and G. Silva-Ortigoza, Class. Quantum Grav., 23, 497-509, 

(2006) . 

[21] C. Kozameh, E.T. Newman and Gilberto Silva-Ortigoza, 
Class.Quantum.Grav.24, 1955-1979, (2007). 

[22] E. T. Newman, Phys. Rev. D, 65, 104005, (2002). 

[23] Ezra T. Newman, Class. Quantum Grav 21 3197-3221,(2004). 

[24] W. Thirring, Principles of Quantum Electrodynamics, (p. 25), Academic 
Press, NY and London, (1958). 

[25] Landau and E. M. Lifschitz, Classical Theory of Fields, Addison- Wesley, 
Reading , Mass., (1962). 

[26] I. Robinson and A. Trautman, Proc. Roy. Soc, London Ser., A 289, 463 
(1962). 

[27] C. J. Talbot, Commun. Math. Phys. 13, 45, (1969). 

[28] C. Kozameh, E.T. Newman and G. Silva Ortigoza, Class. Quantum Grav. 
23, 6599-6620, (2006). 

[29] M. Rivas, J.M. Aguirregabiria, A. Hernandez, Phys. Letters. A, 257, 21-25, 
(1999). 



28 



7 Appendix 



7.1 The Canonical Slicing 

To construct the canonical slicing we begin with the complex UCF 

u = X(t,CX) (104) 
= ^=e°(r) - ^(^(C,C) +£ y (r)l&(C,C) + - 

and write 

t = s + iX (105) 
with s and A real. The cut function can then be rewritten 

u = X{T,CC)=X(s + i\,C,C) (106) 

= xr{s, a, C C) + »X/(s, a, c, 0, 

with real Xi?(s,A,C,C) and x/(s,A, CO- Tne Xfl( s A C,C) and Xi(s, A,C,C) are 
easily calculated from -X(r, £, C) by 

Xfl(«, A, CO = \{X(s + iX, C, C) + X(s + iX,(X)} (107) 

1 



Xi(s, a, c, C) = ^{X(s + iX, c, - X(s + iX, c, c)}. 



By setting 

and solving for 

we obtain the real slicing, 



Xl (s,X,C,O=0 (108) 
A = A(s,C,C) (109) 

u = xr(s,A(s,(,(),(,0. (HO) 



Remark 6. We remark without proof' 61, that using the gauge freedom 
described early for the choice of the parameter r, we can normalize the real 
velocity vector, w^(s) = £fj (s = r), to one, i.e., Vabi'^s^^s) = 1. 



7.2 Lorentzian Tensors 

In Sec. IV we pointed out that certain spin and conformal weighted function on 
the sphere carried finite dimensional representations of the Lorentz group, i.e., 
they carried information about Lorentzian tensor objects. As examples of this 
we will work out two specific cases. 
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Starting with the Lorentz transformation 



ad- h c=l (111) 



c( + d 



e 



(112) 



K = + ^ „ (113) 

[(a( + b) (3C + 6) + « + d) (cC + d)] 

^((■^^(^(CO (114) 
we choose the special transformation 



C = a 2 C (115) 
' a 
; ~~ a 

aa(l + CO 



e iA = ~ = - 
d a 



K = 



[M 2 CC + i] 



Example 1 s = and w = 1 

Applying this special transformation to the invariant subspace of an s = 
and w = 1 quantity, e.g., to u — X{t, £, £), we have 



%>,i)(C,C)= ^7(o,i) (CO (116) 

eulo = kcucq (H7) 
wc,c) = ^(1,^,-^4,^% (us) 

v j 2 v ' l + cc i + CC i + CC 

After using Eq. (|115p in Eq. (|117p and comparing the coefficients of (1,£, C, CC), 
we find that 



C° = i(aa + a" 1 !" 1 )^ + ^(a^ 1 - aa)f (119) 
? = ( l -aa + ^a- 1 )^ 3 + (^"'S" 1 - (120) 



e l -<l 2 = -« 1 -<e a ) (i2i) 



a 



a 

Since = ,can be written as e i(p ,we have a spatial rotation in the (1, 2) plane. 
Then by identifying 

(l-^)- L * =±-(aa + a- 1 a- 1 ) (123) 
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we have the Lorentz transformation 

_ to Hf3 

f° = + ^s^r, (124) 

e 5 - — + — ^V^. (125) 

We see that the special fractional linear transformation £ = a 2 £ corresponds to 
the standard Lorentz transformation with a spatial rotation. 

Example 2 s = — 1 and to = 1 coming from the s = 1 and io = — 3 
isomorphism. 

Applying ( = a 2 £ to the s = — l,w = 1, case, e.g., to the invariant factor 
space of V? = < (1 ,_ 3) = + ^i°« y 2« + - 

we have 

r} { -i,i) = fe(C,C) = e-* x K Vi _ hl) = e-^H(CC) (126) 

™a(C,C) = ^(o,i-C 2 ,-*(i + C 2 ), 2C). 

Comparing the coefficients of (1, £, £ 2 ), we find that 

®,i = \^ + - 4 ( ° 2 - 0_2 ^° a (12?) 



2 



^, 2 = i[a- 2 + o a ]< a + ^[a a - a -Xi 



(128) 



< 2 ^Vii 



If we had identihed with a Maxwell field via 

< = (^ 1 ,i^i°,2,<3) = (E + iB) i , (130) 

then Eq. (|127p would be equivalent to a Lorentz transformation of the Maxwell 
tensor F ab . The six real components of ip^ thus corresponds to a skew-symmetric 
Lorentzian tensor [231. 



7.3 Products of Spin-s Harmonics 

For completeness we give several of the relevant Clebsch- Gordon products that 
were used. We have left out terms with Z-values greater than two. 



31 



Yi kl Y^ = 6 ^li+ l ^5 jl e lke Y? e - U l3 Y^ (131) 
24 „ „ 24, 



Y 2ki Y 2ij - -r-SikSji + -z-SijY 2ik , 



I 2kl x 2lj 



5 „ ? 
£^A V-1^ 12 AV-1 

— g — dikei jf r lf +Y°U Y 2kj> 



yV-i _ ix.. _ '^ r ..,y° _ _LyO 



V 2 V- 1 — ^A V 1 V 1 . y2 yO _ «2\/2 2 

X 2ij x lk ~ ~x°jk x li ~~Q~ eiklI 2jl > " / 2y / lfe — 3 t ikl x 2jli 



V 2 V° - - 24 A V 2 . V 2 V- 1 — 2v ^ ,^ c V 1 4- 6 A V 1 

I 2ij 1 2ml — ~°lj 1 2irm * x 2lm x 2ij — "'"ifflElj/ 1 !/ T y°h x 2mji 



Y^miYii = -Y^Sim — -V2e im fY 2 \f , • Y^Y^ = -Y u 1 Sik + - V2eikfY 2l j , 

yOyO _ ^ 1 o , yOyO _l^A..y° 

7.4 Angular Momentum Ambiguities 

As we mentioned earlier in the text there have been ambiguities, described in 
the literature 0J, in the definition of the asymptotic angular momentum, J . In 
our notation, (omitting the Maxwell field), the ambiguities are in the arbitrary 
choice of the constant p in the expression: 

J* k = -^i>xu\x+P^M^+\^)]k- (132) 

The default choices appear to be either 2, 1, or 0. The present work does not 
influence or help resolve the ambiguity since to second order the expression 

lm[am+^(aa)]k = i^^Vf +f f J >y fc) (133) 

vanishes. 

One might have thought that our flux law, (still omitting the Maxwell field), 
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J k ' = -^^(&v i £+Zfv?)e ljk , (134) 

ji- ^ c \ h o\ 
J --J2G^ 11 

would have an ambiguity in the flux, i.e., a total derivative, arising from the use 
of the chain-rule, e.g., 

This apparent ambiguity disappears since these total derivatives are equivalent 
to the expression in Eq. (|133j) and again vanish identically . 

It is possible that if our calculations were repeated, but done to third order, 
the ambiguities could be resolved. 
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